A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

环境规制强度影响能源效率的门槛效应研究——基于PSTR模型

利用非径向-SBM和CCR模型测算了我国30个省市(暂不包括西藏和港澳台地区)2007-2016年间的环境规制强度以及能源效率,构建面板平滑转移回归模型(PSTR),分析了环境规制强度对我国能源效率的连续非线性变化的影响.研究结果显示,全国环境规制强度对我国能源效率存在非线性效应,二者之间呈倒“U”型关系.环境规制强度小于门槛值0.6376时,环境规制对能源效率产生创新补偿效应,当环境规制强度高于门槛值0.6376时,环境规制对能源效率产生遵循成本效应.而我国东西部地区的环境规制强度与地区能源效率间却存在着与之相反的“U”型效应.     

基于MRS-SJC-Copula模型对A股与港股的动态联动性研究

由于经济增长周期变化,导致不同股票市场存在高低状态转换的现象,在研究不同股票市场之间联动性的研究时,需要考虑股票波动均值和方差两种结构变化。基于具有马尔可夫状态转换的动态SJC-Copula,结合修正ICSS算法对我国内地股票市场和香港股票市场之间的联动性进行方差结构突变点的检验。实证结果表明:国内和香港股票市场之间存在非线性非对称的时变相依性,并持续存在高低两种不同状态的概率转换。股票指数由于动态联动受到负面消息的下跌幅度大于正面消息的变化幅度,且上下尾部均受上期信息的持续影响。“沪港通”、“深港通”、中美贸易战等因素使得其上下尾部发生结构突变,内地与香港股票市场的联动性增大和市场波动幅度趋强。     

Phase transitions in graphs on orientable surfaces

Let be the orientable surface of genus and denote by the class of all graphs on vertex set with edges embeddable on . We prove that the component structure of a graph chosen uniformly at random from features two phase transitions. The first phase transition mirrors the classical phase transition in the Erd?s‐Rényi random graph chosen uniformly at random from all graphs with vertex set and edges. It takes place at , when the giant component emerges. The second phase transition occurs at , when the giant component covers almost all vertices of the graph. This kind of phenomenon is strikingly different from and has only been observed for graphs on surfaces.     

Combining Complex and Radial Slave Boson Fields within the Kotliar–Ruckenstein Representation of Correlated Impurities

The gauge symmetry group of any slave boson representation allows to gauge away the phase of bosonic fields. One benefit of this radial field formulation is the elimination of spurious Bose condensations when saddle-point approximation is performed. Within the Kotliar–Ruckenstein representation, three of the four bosonic fields can be radial while the last one has to remain complex. In this work, the procedure to carry out the functional integration involving constrained fermionic fields, complex bosonic fields, and radial bosonic fields is presented. The correctness of the representation is verified by exactly evaluating the partition function and the Green's function of the Hubbard model in the atomic limit.     

Comparison of the stability of thymine tautomers and the interaction of its tautomers with Na+, K+, Mg2+ and Ca2+ in gas and solvent phases

The present study compared the interactions among Na ⁺, K ⁺, Mg²⁺ and Ca²⁺, thymine and its tautomers in the gas and solvent phase, an interaction dependent upon the electronic construction of the tautomers. Three types of cation interaction with thymine and its tautomers were observed. In the first one, the metal cations interacted with a lone pair of nitrogen or oxygen of the tautomers. In the second type, there was an interaction among the cations, nitrogen and oxygen at the same time; the last one was that of cations with the electron density of thymine π-system, where the cations were perpendicular to the ring of thymine. The interaction of metals cation with tautomers was studied in the gas and solvent phases; a comparison was then made between interactions in two phases. The interaction energy for all complexes indicated the stability of complexes, an energy which was higher in Ca²⁺ and Mg²⁺ compared with Na⁺ and K⁺. Concerning K⁺ and Na⁺, the stability of all complexes of tautomers was greater than that of thymine complexes; however, the stability of certain Ca²⁺ and Mg²⁺ complexes was lower than the complexes of thymine.     

Quantum to classical crossover under dephasing effects in a two-dimensional percolation model

Scaling theory predicts complete localization in d = 2 in quantum systems belonging to the orthogonal class(i.e., with timereversal symmetry and spin-rotation symmetry). The conductance g behaves as g^exp(-L/l) with system size L and localization length l in the strong disorder limit. However, classical systems can always have metallic states in which Ohm’s law shows a constant g in d=2. We study a two-dimensional quantum percolation model by controlling dephasing effects. The numerical investigation of g aims at simulating a quantum-to-classical percolation evolution. An unexpected metallic phase, where g increases with L, generates immense interest before the system becomes completely classical. Furthermore, the analysis of the scaling plot of g indicates a metal-insulator crossover.     

Tolerance factor and phase stability of the garnet structure. Corrigendum

An error in an equation in the paper by Song et al. [ Acta Cryst. (2019), C 75 , 1353–1358 ] is corrected.     

Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.